Cluster Percolation and Chiral Phase Transition 
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The Meron Cluster algorithm solves the sign problem in a class of interacting fermion lattice 
models with a chiral phase transition. Within this framework, we study the geometrical features of 
the clusters built by the algorithm, that suggest the occurrence of a generalized percolating phase 
transition at the chiral critical temperature in close analogy with Fortuin-Kasteleyn percolation in 
spin models. 

PACS numbers: 71.10.Fd, 64.60.Ak 
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A fundamental difficulty in the Monte Carlo study of 
fermion lattice models is the sign problem due to the 
fluctuating sign of the statistical weight of fermion con- 
figurations fy] . Recently, the Meron Cluster algorithm 
(MCA) (|||has been proposed as an effective solution 
to the sign problem in a class of interacting models. In 
particular, here, we focus on a 2 + 1 dimensional model 
with a second order phase transition associated to the 
dynamical breaking of a discrete chiral symmetry [Q . 

Like all cluster algorithms for lattice models, MCA de- 
fines clusters of sites used as effective non-local degrees of 
freedom to update configurations without critical slow- 
ing down. For a given observable, the sign problem is 
cured by restricting the Monte Carlo sampling to specific 
topological sectors that give contributions not canceling 
in pairs due to the fermion sign. The relevant topolog- 
ical charge is the so-called meron number that we shall 
define later. The rules that assure convergence to the 
correct Boltzmann equilibrium distribution determine a 
well defined cluster dynamics. From the study of lat- 
tice spin models we know that this artificial dynamics 
can be surprisingly rich; in fact, experience in that field 
suggests the existence of a purely geometrical phase tran- 
sition concerning the algorithm clusters and underlying 
the physical thermal transition 

The simplest example of such a scenario is the 2D Ising 
model where clusters can be defined in a natural way 
as sets of nearest neighboring aligned spins and admit 
a physical interpretation as real ordered domains. At 
the thermal transition temperature these clusters per- 
colate ||], but since the critical exponents do not coin- 
cide with the thermal ones J^|, a complete equivalence 
of the two transitions cannot be claimed. In three di- 
mensions, the comparison is even worse and also the two 
critical temperatures are slightly different pc| . To find 
a geometrical transition occurring at the thermal criti- 
cal point with the same critical exponents, generalized 
clusters must be defined jn]|, like the Fortuin-Kasteleyn 
bond clusters and their extensions |^|| . The equivalence 
between the thermal phase transition and a suitable per- 
colative process can then be extended to models with 
continuous rotational or gauge invariance |jj . 



Similar investigations lack for fermionic models with 
sign problems because cluster algorithms have not been 
available until MCA. Here, for the first time, we aim 
to the identification of a geometrical transition in the 
MCA dynamics and look for critical phenomena defined 
in terms of cluster shapes. On the other hand, we must 
keep into account at least the global configuration signs 
because they are the only memory of the fact that the 
model is fermionic and allow to tell it from its bosonized 
counterpart free of sign problems. Besides, apart from 
the sign problem, the rules to build clusters in MCA are 
not precisely the same as for Fortuin-Kasteleyn clusters 
and the existence of a transition is non trivial. 

The model we study describes relativistic staggered 
fcrmions, hopping on a 2 + 1 dimensional lattice with 
L 2 spatial sites, described in B|. The Hamiltonian is 
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where n x = c x c x is the occupation number at site x, r\ x .i 
are the Kawamoto-Smit phases r\ Xi \ = 1, t\ x ,i — i~^) Xl 
and the operators c, c* obey standard anti-commutation 
relations {c x ,c y } = 0, {c x ,c y } = 8 x _ y . In the follow- 
ing we shall consider the case G = 1 and adopt periodic 
boundary conditions. 

The partition function Tr e~@ H can be computed by 
Trotter splitting that maps the quantum model to a sta- 
tistical system on a 2 + 1 dimensional lattice with L 2 xT 
sites. The temporal lattice spacing is e = 4/3/T. The 
limit e — > must be taken at fixed (3 with T — > oo. In 
practice, we shall present results obtained at the fixed 
value T = 40 where (3 can cover the transition point with 
reasonably small e 0. 

Each configuration is specified by the occupation num- 
bers n = {n x t } and carries a sign cr(n) = ±1, source of 
the sign-problem. To update a configuration, sites are 
clustered according to definite rules depending on f3 and 
discussed in details in ||. Each cluster is then indepen- 
dently flipped: with probability 1/2 we apply the global 
transformation n x t — > 1 — n x t to all of its sites. Clusters 
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whose flip changes er(n) are defined merons. 

The chiral phase transition can be analyzed by study- 
ing the asymptotic volume dependence of the susceptibil- 
ity x- Defining the chiral condensate in the configuration 
n = {n X}t } as 



x.t 
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the chiral susceptibility \ is given by 
1 {{Z(nf) q(n)) 
X /3i 2 <a(n)) ■ 1 > 

An improved estimator of x f ree °f s i§ n problems can 
be built by writing Z(n) as a sum over clusters Z(n) — 
Zq and taking its average over 2^° possible flips, 
where Nc is the number of the clusters. Then, \ gets 
contributions from sectors with mcron number N = 0, 2: 
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where, for N = 2, G\, C% are the two merons. 

From numerical simulations, apart from rather small 
scaling violations, x obeys the Finite Size Scaling (FSS) 
law x = L 1 f x {L((3 — (3 t h)) with the exponent 7 = 7/4, 
characteristic of the 2D Ising universality class Jl2| . 
Numerical simulations locate the chiral (3 at (3th = 
2.43(1) |D (the subscript "ift" stands for HhermaT). 

To detect a possible purely geometrical transition, we 
study quantities Q that depend only on the cluster shape 
and not on their internal occupation numbers. Since clus- 
ter flips do not change Q and meron flips change the sign 
of <r(n), the improved estimator of Q is 
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that is the average restricted to the zero meron sector. 

The simplest set of geometrical quantities that we can 
study are the moments M n = E {J2c \C\ n ) °f the nor ~ 
malized cluster size \C\ — Vol(C)/T where Vol(C) is the 
number of sites in C. 



FIG. 1. FSS analysis of M„ and E (max |C|). 

We perform simulations to compute numerically 
{M n } n= 2,3A and also E (max \C\). We work on lattices 
with L = 8, 10, 12, 14, 16 in the range 1.5 < /3 < 3.0 and 
presented about 1.5 • 10 5 measures per point. Fig. (|l|) 
shows the numerical data supporting as a first result the 
remarkable validity of the empirical scaling relations: 



(3 n / 2 M n = L np f n {y), 
(3 1/2 E (max|C|) = If h(y) 



n > 2, 
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in terms of the scaling variable y = L{j3 — (3 C ). This 
scaling behavior defines an order parameter of the geo- 
metrical phase transition. The results for the exponent 
p and the critical (3 C are p = 1.71(5) and (3 C = 2.42(5). 
Within errors, p is consistent with the exact 2D Ising 
exact value 7 = 7/4. We also find (3 C ~ (3th with an 
accuracy that can be appreciated by looking at Fig. (||). 
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FIG. 2. Data crossing for Kh with different L. By virtue 
of the scaling law Eq. at (3 = (3 C the quantity /3M 2 L~ 2p 
is independent on L. 

The ratios R n = M 2 1/2 /Jl4 /n for n = 3,4 and R x = 
,1/2 
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Ma'VEoCmax |C|) are shown in Fig. (|). 
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FIG. 3. Ratios R n and Roc- Mn is well approximated 
by the contribution of the largest cluster beyond the critical 
point. 

The plots strongly indicate that R n are independent on 
L at fixed y in agreement with Eq. (^). 

It can be checked that, at j3 = (3 C , M n receives the main 
contribution from a small set of large clusters growing like 
L p , and the fact that the ratios R n (y) — ► 1 as y grows, 
simply means that this contribution is more and more 
dominant. We briefly summarize this behavior by saying 
that the clusters are percolating. Following |l3[ , further 
information on the critical ensemble can be obtained by 
studying the cluster distribution 

n s = Eo(# number of clusters C with Vol(C) = s). (8) 
In Fig. (|]) we show that the simple law 

n, = L->f (±) , (9) 

is well satisfied for s/L p ^> 1 . 
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FIG. 4. FSS plot of n s . The curves are obtained after av- 
eraging over blocks of 10 subsequent sizes. In the region with 
L p < s < 40 L p , n 3 decreases roughly algebraically; beyond 
the knee at s ~ 40 L p the distribution falls down quickly. The 
critical region is dominated by the clusters at the right edge 
of the plot. 

The critical cluster distribution n s decreases not faster 
than algebraically with s until s ~ 40 L p where the final 
large cluster tail is reached and n s falls down quickly. 

If one takes into account that the leading contribu- 
tion to {M n }„>2 actually comes from the region where 
s/L p ^> 1, observing that M n = J2 s snris aim using 
Eq. (^) we obtain a L-dependence consistent with Eq. (j^). 

There are many small clusters and one can check that 
the normalized number of clusters (3 ~Eq(Nc) / L 2 depends 
mildly on L for a wide range of (3 and L. 

It is interesting to analyze what happens to the typi- 
cal cluster configurations when (3 is gradually increased 



toward the critical point. At small j3, almost all bonds 
that build the clusters are set in the temporal direction. 
In physical terms the fermions hop from the initial po- 
sitions to neighboring ones with small probability. Sites 
simply tend to cluster in L 2 straight vertical lines with 
T sites. As /3 — ► /3 C , the vertical clusters start merging 
and breaking and form complicated structures with large 
dispersion in the cluster size distribution. This process 
can be seen in Fig. (||) : 
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FIG. 5. Two elementary clusters with equal dimension T, 

first merge into a single cluster, then break apart in two clus- 
ters with very different dimension. 

two vertical clusters with Vol(C) = T undergo a two step 
process allowed by the cluster rules 0. In the end, they 
give rise to a large cluster with Vol(C) = 2T — 4 and 
a small one with 4 sites. After many processes of this 
kind we find a few large clusters surrounded by a gas of 
smaller clusters. In Fig. ( |^) 




FIG. 6. A typical largest link configuration for the 8 2 x 40 
lattice at /3 = j3 c . 

we show a typical largest cluster obtained at j3 = 2.5 on 
the 8 2 x 40 lattice. 

Geometrical quantities that can measure this disper- 
sion effect are the cumulants {G n } n >i of the cluster size 
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distribution defined as 
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In our study we do not consider the cumulant G±, 
whose behavior is determined by the the contributions 
from small clusters, with s -C L p . 

The next cumulant is the variance 
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Our data support a FSS law of the form 
f3G 2 = Lp'g 2 (y), 



(11) 
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with p 1 — 1.51(4), which is compatible with 27 — 2. This 
exponent can be explained taking account that, analo- 
gously to the moment M2, the contributions from small 
clusters are negligible and that, within errors, \C\ 2 ~ 
L 21 and N c ~ L 2 . 

As we remarked above, when (3 — ► 0, the clusters are 
simple vertical lines and G2 — > like x does. This fact 
suggests that it could be interesting to study a possi- 
ble relationship between these quantities at different (i 
values. We remark that, in principle, G2 and x have dif- 
ferent topological origins: x is calculated on zero and two 
merons sector, while G 2 only on zero-meron sector. Nev- 
ertheless, we find that the following empirical relation 



X *l 0.15(1) • (fi G 2 )~- 
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holds for a wide range of parameters (L, (3) as shown in 
Fig. ©. 
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FIG. 7. Comparison between \ ( em Pty symbols) and 
0.15 • (/3G 2 ) 7/6 (full symbols). 

In practice, in the region that we have explored, the ratio 
x/(/3G2) 7 / 6 is consistent, within errors, with a constant. 
This result signals an unexpected correlation between dif- 
ferent topological sectors. In particular, Eq. ( |T3| ) could 



be relevant in the construction of a purely geometrical 
definition of x- A FSS study of the higher cumulants G3 
and G4 shows similar scaling laws, but with exponents 
that are not in simple relation with 7 and should in prin- 
ciple be matched to the anomalous dimensions of higher 
operators in the 2D Ising universality class. 

To conclude, we have examined the critical behavior of 
the clusters that arise in the application of the Meron al- 
gorithm to a fermion model in 2 + 1 dimensions. We have 
found simple FSS laws that we have explained in terms of 
a percolative process occurring at the chiral critical tem- 
perature. Our data support the results p ~ 7, p' ~ 27—2, 
as well as the empirical relation Eq. (|l3| ) that shows a 
close correlation between physical and geometrical quan- 
tities. 
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